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^H | Abstract: In this paper, we show how a power-law correction to the Einstcin- 

Hilbert action provides a viable modified theory of gravity, passing the Solar- 
System tests, when the exponent is between the values 2 and 3. Then, we imple- 
ment this paradigm on a cosmological setting outlining how the main phases of 
the Universe thermal history are properly reproduced. 

As a result, we find two distinct constraints on the characteristic length scale 
of the model, i.e., a lower bound from the Solar-System test and an upper one 
0^ i by guaranteeing the matter dominated Universe evolution. 
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1 Basic statements 

From the very beginning, the possibility to riformulate General Relativity by using a generic 
function of the Ricci scalar (see, for example, [1 for a recent review and references therein) 
appeared as a natural issue offered by the fundamental principles established by Einstein. 
However, it is important to remark that any modification of the Einstein-Hilbert (EH) La- 
grangian is reflected onto a deformed gravitational-field dynamics at any length scale investi- 
gated or observed. Thus, the success of such f(R) gravity in the solution of a specific problem 
has to match consistency with observation in different length scales [31 H] . A viable self- 
consistent model can be often obtained at the price to consider a generalized gravitational 
lagrangian containing a large number of free parameters. Nevertheless, the wide spectrum of 
possible choices for f(R) can appear as a weakness point in view of the predictivity of the 
theory, because a significant degree of degeneracy is expected in the model. 

Here, we consider an opposite point of view, by studying the viability of a power-law 
correction to the EH action having a single free parameter (a length scale) once the power- 
law exponent is fixed. We investigate the implementation of the Solar-System test to our 



model [5| and then we pursue a cosmological study of the resulting modified Fricdmann- 
Lemaitre-Robertson- Walker (FLRW) dynamics. As expected, this scenario gives us a rather 
stringent range of variation for the free length scale where searching for new gravitational 
physics. 

2 Non-analytical power-law f(R) model 

In this paper, we consider the following modified gravitational action in the so-called Jordan 
frame 

S=-^Sd i x^—gf{R), f(R) = R + qR n , (1) 

where n is a non-integer dimensionless parameter and q < has dimensions of [X] 2 ™ -2 (in 
the equation above x = 87rG, using c = 1 and G being the Newton constant, moreover, the 
signature is set as [+,—,—, — ]). Such a form of f(R) gives the following constraints for n: 
if R > 0, all rt-values are allowed; if R < 0, the condition n — £/(2m + 1) must hold (where, 
here and in the following, m and £ denote positive integer). It is straightforward to verify 
that S in eq.tUJ) is non-analytical in R = for non- integer, rational n, i.e., it docs not admit 
Taylor expansion near R = 0. 

Let us now define the characteristic length scale of our model as 

L q (n) = | g |V( 2 «- 2 ) , (2) 

while variations of the total action Stot — S + Sm (where Sm denotes the matter term) with 
respect to the metric give, after manipulations and modulo surface terms: 

f'R»,v - \ 9ixvf - V^V^ /' + g MV a /' = x T iiu , (3) 

where T^ u is the Energy- Momentum Tensor (EMT). Here and in the following (...)' indicates 
the derivative with respect to R, □ = g pa V P V a and V M or (...) ; denotes the covariant 
derivative (Greek indices run form to 3). 

We can gain further information on the value of n by analyzing the conditions that allow for 
a consistent weak- field stationary limit. Having in mind to investigate the weak field limit of 
our theory to obtain predictions at Solar-System scales, we can decompose the corresponding 
metric as <? M „ = n^ + h^ v , where /i M „ is a small (for our case, static) perturbation of the 
Minkowskian metric rj pl/ . In this limit, the vacuum Einstein equations read 

JV - h^ R - nqiR 11 - 1 )-^ + nqri^nR 71 - 1 = , R = inqnRJ 1 - 1 . (5) 

The structure of such field equations leads us to focus our attention on the restricted region 
of the parameter space 2 < n < 3. This choice is enforced by the fulfillment of the conditions 
by which all other terms are negligible with respect to the linear and the lowest-order non- 
Einsteinian ones. 

3 Viability of the theory: the Solar-System test 

From the analysis of the weak-field limit in the Jordan frame, i.e., eqs.([5]), we learn the 
possibility to find a post-Newtonian solution by solving eqs.([5]) up to the next-to- leading 
order in h, i.e., up to 0(h n ~ l ), and neglecting the 0(h 2 ) contribution only for the cases 
2 < ii < 3. These considerations motivate the choice we claimed above concerning the 
restriction of the parameter n. 

The most general spherically-symmetric line element in the weak-field limit is 

ds 2 = (1 + ®)dt 2 - (1 - ^)dr 2 - r 2 dn 2 , (6) 



where $ and \fr are the two generalized gravitational potentials and dVL 2 is the solid angle 
element. Within this framework, the modified Einstein equations ([5]) rewrite 
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where (...). denotes ordinary differentiation. The system above is solved by 
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where the integration constant 5 has the dimensions of [L] and the dimensionless integration 
constant a can be set equal to zero without loss of generality. The integration constant A has 
the dimensions of [_L]( 2n_2 )/( 2_n ), and <&„ and \&„ are dimensionless, accordingly. Moreover, 
one can check that $„ and ^> n are well-defined only in the case n = (2m + l)/£ while we 
get A > since we assume q < 0. In agreement to the geodesic motion as expanded in the 
weak field limit, the integration constant 6 results equal to 6 = — r s , where r s = 2GM is the 
Schwarzschild radius of a central object of mass M. 

The most suitable arena to evaluate the reliability and the validity range of the weak-field 
solution © is the Solar System [2HJ. To this end, we can specify eas. (|8b)l - ([5c)) for the typical 
length scales involved in the problem and we split $ and \& into two terms, the Newtonian 
part and a modification, i.e., 



* = $» + *M 



-rg/r + ^r/L,) 2 ^/^- 2 ), 
-rg/r + vMr/L^-DA"- 2 ), 



(9a) 
(9b) 



here, the integration constant 5 of eas. ((8b|> - (|5c)l is 5 = — r = 2GM (M being the Solar 
mass). While the weak-field approximation of the Schwarzschild metric is valid within the 
range rg <C r < oo because it is asymptotically flat, the modification terms have the peculiar 
feature to diverge for r — > oo. It is therefore necessary to establish a validity range, i.e., 
fMin <C r -C rji./ai, related to n and L q , where this solution is physically predictive [BJ. 

Since we aim to provide a physical picture at least of the planetary region of the Solar 
System, we are led to require that <&m and ^ m remain small perturbations with respect to 
$at and *&n, so that it is easy to recognize the absence of a minimal radius except for the 
condition r 3> rg. The typical distance L Q corresponds to the request 



\*n(L*)\ ~ \$ M (L*)\ 



|*iv(£ )| ~ I*m(£ )| 



(10) 



For rg -C r <C L Q , the system obeys thus Newtonian physics and experiences the post- 
Newtonian term as a correction. Another maximum distance Lq* can be defined, according 
to the request that the weak-field expansion ([6J should hold, regardless to the ratios <&m/®n 
and ^m/^n- L q * results to be defined by 



|cML *)| « \$ M (L%)\ ~ 1 



1*^(^)1 « |*m(£ *)| ~ 1 



(11) 



We remark that Lq and Lq are defined as functions of n and L qi i.e., 

n—2 2 yi — 2 n — 2 

^ ~ |r|/$ n |^=4L|™- 4 , L- ~ Vl^l 5 ^ > ( 12 ) 

and it is important to underline that, for the validity of our scheme, the condition Lq 3> Tq 
must hold, i.e., L q > rg|$„|(™- 2 )/( 2 "- 2 ). 

Neglecting the lower-order effects concerning the eccentricity of the planetary orbit, we 
can deal with the simple model of a planet moving on circular orbit around the Sun with an 
orbital period T given by T — 2n(r/a) 1 ' 2 (a = (d<&/dr)/2 being the centripetal acceleration). 
For our model, from eas. (|8bp . we get 
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We now can compare the correction to the Keplerian period Tk = 2ttt 3 /' 2 (GMq)~ 1 / 2 , with 
the experimental data of the period T exp and its uncertainty 5T exp . We then impose the 
correction to be smaller than the experimental uncertainty, i.e., 



Olexp > | J K — -Ln 
■*■ exp *- j( 



I $ »I^(^ 2 )/(^V- 2 ), (14) 



where rp is the mean orbital distance of a given planet from the Sun. 

Let us now specify our analysis for the example of the Earth [5] . In this particular case, 
T exp ~ 365.2563 days and ST exp ~ 5.0 • 1CT 10 days (with r P ~ 4.8482 x 10~ 6 pc). This way, 
for the Earth, we can get a lower bound L q > L^ n for the characteristic length scale of our 
model, as function of n, i.e., 



I*J n-1 



3w-4, _™ 



2n-2 



i7(») = 1-3689 x 10" ^ r£-* '-' , (15) 

H L rg n — 2 J 

where $„ is defined in eq. (|8b[) and L^ m ~4x 10~ 3 pc, for a typical value n ~ 2.66. We 
remark that L^ n , by virtue of eq. (|8b|) . is defined only for n = (2m + \)ji. 

Our analysis clarifies how the predictions of the corresponding equations for the weak-field 
limit appear viable in view of the constraints arising from the Solar-System physics. Indeed, 
the lower bound for L q does not represent a serious shortcoming of the model, as we are 
going to discuss in Sec[6j where a plot of L q Q n {n) and of Lg(£ g ) and Lq (L q ) will be also 
addressed. 

4 Cosmological implementation of the f(R) model 

In order to study how our f(R) model affects the cosmological evolution, we start from the 
modified gravitational action ([l) and we assume the standard Robertson- Walker (RW) line 
element in the synchronous reference system, i.e., 

ds 2 = dt 2 - a(tf [ dr 2 /(l -Kr 2 )+ r 2 dil 2 ] , (16) 

where a(t) is the scale factor and K the spatial curvature constant. Using such expression, the 
00-component of eq.Q results, for symmetry using the Bianchi identity, the only independent 
one and it writes as 

f'Roo-hf + 3(d/a)f"R = X Too. (17) 



where the dot indicates the time derivative. We assume as matter source a perfect-fluid EMT, 
i.e., T^ v = (p + p) u^u v —pg^ v , in a comoving reference system (thus Too = p), where p is the 
thermostatic pressure, p the energy density and u^ denotes the 4-velocity. The 0-component 
of the conservation law, i.e., T" v — with v — 0, assuming the equation of state (EoS) 
p = wp, gives the following expression for the energy density: p = p a [a / 'a ]~ 3 ( 1+K '). 
Using now / = R + qR n with q < 0, we are able to explicitly write ea. (fT7|) : 

2xa x - 3w + 6 n nqa 5 -' 2n a(-K-a 2 -aa) n - 1 + 

+ a 2 [-6K-6a 2 + 6 n q a 2 ^~ n ^ (-K - a 2 - a a) n ] + (18) 

+ 6"(ii-l)n(jda 2(2 "" ) (-K - a 2 - ah) n - 2 [- 2d 3 - 2Kb, + aaa + a 2 a] =0 , 

where x = XPo a o ■ Let us now assume a power-law a = a [t/t ] x for the scale factor and, 

for the sake of simplicity, we set a — a t~ x (clearly, [a] — [L 4 ~ x ]). Here and in the following, 
we use the subscript (...) to denote quantities measured today. In this case, eci.(fT8f can be 
recast in the form 

-6a 2 Kt 2x -6a i x 2 t 4x - 2 + qa i t 4x (c 1 t- 2 ~6a- 2 Kt- 2x y+ 2 X a 4 - 3w t^ 3 ^ = 

where d = 6x(l - 2x), C 2 = (x(2n - 1) - 1), C 3 = xa 2 (x + 2n- 3){2x - 1), and 
C A = xa 2 {2x~ 1). 

4.1 Radiation-dominated Universe 

Here, we assume the radiation-dominated Universe EoS w — 1/3 (p ~ a -4 ). In the following, 
we will discuss the three distinct regimes, in the asymptotic limit as t — > 0, for x < 1, x > 1 
and x = 1, separately. 

In the case x < 1, all terms containing explicitly the curvature K of ea. dT^ll results to be 
negligible for t — > and asymptotic solutions are allowed if and only if x sj n/2 which, in 
the case 2 < n < 3 we are considering, is always satisfied. The leading-order term of ea. fTSj) 
writes as 

qa 4 C? [ 1 - (C 3 /Cl)nx 2 a 2 } t 4x ~ 2n = , (20) 

andx = 1/2 and a; = [2 + 3n-2n 2 ± (4 + 8n + n 2 - 12n 3 + 4n 4 ) 1 / 2 ]/2n are the solutions. Such 
second expression results to be negative or imaginary for 2 < n < 3 and must be excluded. 
Thus, the only solution for x < 1, in the asymptotic limit for t — > 0, is the well-known 
radiation dominated behavior a ~ t 1 ' 2 . In the other two cases, i.e., for x ^ 1, it is easy to 
recognize that no asymptotic solutions are allowed. Therefore, the approach to the initial 
singularity is not characterized by power-law inflation behavior when spatial curvature is 
non- vanishing. 

Let us now assume a vanishing spatial curvature in ea. (|19[) . In can be show how, for 
K = 0, the radiation-dominated solution with w — 1/3 and x = 1/2 is an exact solution 
(non-asymptotic and allowed for all n-values) giving p = 3/(4%^), matching the standard 
FLRW case. In the case x > 1, the leading-order terms of eq. tfTO]) read, for t — > and K = 0, 



qa A C? [ 1 - {C 3 /Cl)nx 2 a 2 } t 4x ~ 2n + 2 X = . (21) 

Three distinct regimes have to be now separately discussed. For x > n/2, the leading order of 
the equation above does not admit solutions since it writes simply 2\ — and, for x < n/2, 



the solutions of ea. (j2Tj) are those obtained in the case for x < 1. Instead, for x — n/2, and 
defining H a — (n/2)/t , one gets 

P° = I TW 1 ' ~P^ = T C 1 ~ n) ( " _1) n n (n(4 + (6 - 5n)n) - 4)(n/2) 2 - 2 ", (22) 

where we have introduced the dimensionless parameter q = H 2n ~ 2 q. We remark that the 
constraint n = (2m + l)/£ (which is in agreement with respect to the one obtained from Solar- 
System test) must hold in order to have p > since we have assumed q < and therefore 
q < 0. The function p results to increase as n goes from 2 to 3 and, in particular, one can 
get 216 < p < 21 024. Finally, for x = 1 and K = 0, eq.ip)) reads [1 - n{2n - 2)] £ 4 ~ 2n = 0, 
giving n = [1 ± v3]/2. As the previous case, the regime x = 1 does not admit solutions in 
the region 2 < n < 3. 

4.2 Matter-dominated Universe 

Let us now study the matter-dominated Universe EoS w = (p ~ a -3 ). As previously done, 
we analyze the three distinct regimes for x < 1, x > 1 and x = 1, and, in the limit for 
i — > oo, it is easy to recognize that there are no power-law solutions in all these cases for 
K 7^ 0. Setting K = 0, the x ~^ 1 regimes do not provide any power-law form for cosmological 
dynamics either. On the other hand, for x < 1 and assuming zero spatial curvature in ea. ([T9| . 
we get the following equation: 

[-6x 2 a 4 ] t 4x - 2 + 2x5 t x = a 4 qC?[-l + C 3 U 2 nx 2 /Cl] t ix - 2n . (23) 

Since 4x — 2 > 4x — 2n, the term on the right hand side can be neglected in the limit of large t 
and the equation above admits three distinct situations: x < 2/3, x > 2/3 and x = 2/3. Both 
cases with x ^ 2/3 do not admit solution. The case x = 2/3 admits instead an asymptotic 
solution for t — > oo. In fact, ea.([23| reduces to 8a 3 = 6\ and the FLWR matter-dominated 
power-law solution a — at 2 / 3 is reached setting p = 4/(3x* 2 )- 

In conclusion, we can infer that, for f(R) = R + qR n , the standard matter-dominated 
FLRW behavior of the scale factor a ~ t 2 / 3 is the only asymptotic (as t — > oo) power-law 
solution. 

4.2.1 Range of t-values: 

As shown above, the matter dominated solution a ~ £ 2 / 3 is obtained for K = and 
asimptotically as t — > oo. In order to neglect all the fiT-terms in our f(R) model, we start 
directly from the expression of the Ricci scalar [7]. Using a power-law scale factor, we get 
the t-range (if x ^ 1/2 and x < 1) 

t « | [x(2x - l)]/[K/a 2 } \ 1/{2 - 2x) . (24) 

For the matter-dominated era and using standard cosmological parameters [5], one can get 
the upper limit K/a 2 < 0.006 (H ) 2 / 3 , to estimate the value of K/a 2 . Thus, setting x = 2/3, 
we get the bound t <C 235/ H , independently of the form of f(R). 

At the same time, if we set x — 2/3, the asymptotic solution p — 4/(3%£ 2 ) is reached 
neglecting the right hand side (-C 1) of eg. (f2"5)l . i.e, if t is constrained by the following lower 
limit (we remind that q — H 2n ~ 2 q) 

£» n(n,q )/H , p{n,q a ) - | q [ - (4/3)" + 2( 2n+1 )3^ n(2n - 7/3)] | 1/2(n_1) . (25) 



Let us now recall that the matter-dominated era began, assuming H~ 4 ~ 4.3 x 10 s, at 
tEq — 5.1 x 10~ 6 /if Q . In this sense, we can safely assume p(n, q ) ^ 5.1 x 10 -8 , which implies 
an upper limit for |g |, i.e., \q \ < \q \ Max , where 

\q \ Max {n) = [5.1xl0- 8 ] 2(1_n) |-(4/3)" + 2( 2n+1 )3- n n(2n-7/3)f 1 . (26) 

It is easy to check that the function \q \ Max (n) is decreasing as n goes from 2 to 3, in 
particular, one gets: 10 -16 < \q \ Max < 10~ 31 . 

5 The inflationary paradigm 

After discussing the power-law evolution of the Universe proper of the radiation- and matter- 
dominated eras, we now analyze the inflationary behavior characterizing the very early dy- 
namics (for an interesting approach to the inflationary scenario within the modified gravity 
scheme, see [SIHD])- In this respect, we hypothesize an exponential behavior for the scale 
factor of the Universe a = a e s (' _ *°) = a e st , where s > and a = a e~ st °. In the following, 
we concentrate the attention on the solution for vanishing spatial curvature K = and, in 
this case, ea. ([T51) rewrites as 

a 4 e 4st [q(~12) n s 2n (l - n/2) - 6s 2 ] + 2x{ae st ) 1 ~ 3w = . (27) 

Let us now assume w = —1 (i.e., p = pi = const.) during inflation. Using the definition 
q = H 2n ~ 2 q, the equation above reduces to 

[(-1)"12" q (l - n/2)] s 2n - 6s 2 + k = , (28) 

where k = 2\piH~ 2 and s is a dimensionless parameter defined as s = s/H . Since H de- 
notes the Hubble parameter measured today and estimating H] — yXPifii {i.e., accordingly 
to its Friedmannian value) during inflation [7], one can obtain k ~ H 2 / H 2 ~ 10 100 . For such 
values, it is easy to realize that considering the case n — 2£/(2m + 1), the equation above 
does not admit real solution, thus we now discuss, consistently with the previous analyses, 
only n = (2m + l)/{2t + 1). 

In order to integrate ea. (|28p . we focus on a particular value of the power-law f(R) exponent, 
e.g, n = 29/13 ~ 2.23. Using eq. p^)) . for this value of n one obtains that it can be safely 
considered L^ n ~ 1.44 x 10" 5 pc and, having in mind that L q — Iqol 1 ^ 2 "^ 2 " 1 /H with 
H- 1 ~ 4.2 x 10 9 pc, we get \q \ > \q„\ Mm - 2.56 x 10- 36 . 

Let us now fix the parameter q to a reasonable value like q* ~ — 10 3 |g o | " (such as- 
sumption will be physically motivated in the next Section). In this case, the solution of 
eq . (l28|) is s ~ 2.45 x 10 29 . This analysis demonstrates that an exponential early expansion 
of the Universe is still associated to a vacuum constant energy, even for the modified Fricd- 
mann dynamics. However, we see that the rate of expansion is significantly lower than the 
Friedmann-like one of about a factor in s of 10 20 . Although our estimation relies on the 
Friedmannian relation between Hi and pi (the latter is taken of the order of the Grand Uni- 
fication energy-scale) , nevertheless the values of s remains many order of magnitude below 
the standard value ~ 10 50 even if we change Hj for several order of magnitude. Despite 
this difference, it is still possible to arrange the cosmological parameter in order to have a 
satisfactory inflationary scenario, as far as we require a longer duration of the de Sitter phase. 

6 Physical remarks 



As already discussed in Sec. 2, the parameter q has dimension [L] 2n 2 . We have therefore 



defined a characteristic length scale of the model as L q (n) — Igj 1 /' 2 ™ 2 \ Assuming f(R) 



corrections to be smaller than the experimental uncertainty of the orbital period of the Earth 
around the Sun, the lower bound (|15[) for L q {n) was found. In order to identify the allowed 
scales for our model and in view of the upper constraint on the parameter q a 
derived in the cosmological framework, we can now define the upper limit for L q (n) as 



H 2 n ~ 



j Max 



(n) 



J*. 



Maxll/(2n- 



2) /H Q , 



(29) 



which, considering eo. (l2"6"l) . yields to the constraints 65.59 pc < L q ax < 78.37 pc, for 
2 < n < 3. Assuming H' 1 ~ 4.2 x 10 9 pc, the two bounds for the characteristic length 
scales here discussed, i.e, eq. (fl~5|) and cq. ((29|) . are plotted in FigfljA). At the same time 





Figure 1: Panel A: 



t Mir 



of eq.(fT5 



and Lf ax 



of eq. (f2l?|) . The gray zone represents the 



allowed characteristic-length scales of the model. We stress that L^q™ is defined 
only if n — [2m + !)/(, as represented by the dotted line. Panel B: L & and L Q * 
of eq. ([j"2j) . The gray zone represents here the allowed distances for the model. 



two other typical lengths have been outlined in eq. (fT2")) for the Solar System. L Q represents 
the minimum distance to have post-Newtonian and Newtonian terms of the same order. 
While Lq* was defined according to the request that the weak-field expansion holds. Setting 
now n = 23/9 ~ 2.55, one can show from eQ. flT5"j) and eq. (j2T)]) that the allowed scales are 
0.0013 pc < L q < 71.72 pc. In this range, L Q and Lq can be plotted as in Fig[HB). 

Summarizing, our analysis states a precise range of validity for the power-law f(R) model 
we consider. Indeed, for a generic value of n (i.e., not close to 2 or 3) the fundamental length 
of the model is constrained to range from the super Solar-System scale up to a sub-galactic 
one. Therefore, in agreement to eo. (|9ap . we have to search significant modification for the 
Newton law in gravitational system lying in this interval of length scales, like for instance, 
stellar clusters. 
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